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^ Abstract 
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We use Monte Carlo simulations to study pure 2D Euclidean quan- 
Q,^ tum gravity with i?^-interaction on spherical topologies, employing 

Regge's formulation. We attempt to measure the string susceptibility 
exponent 7str by using a finite-size scaling Ansatz in the expectation 
0^ value of R^, as has been done in a previous study by Bock and Vink 

^ (Nucl. Phys. B438 (1995) 320). By considerably extending the range 

and statistics of their study we find that this Ansatz is plagued by large 
systematic errors. The specific string susceptibility exponent 7^^^. is 
found to agree with theoretical predictions, but its determination also 
is subject to large systematic errors and the presence of finite-size scal- 
^ ing corrections. To circumvent this obstacle we suggest a new scaling 

^ ^ Ansatz which in principle should be able to predict both, 7str and 7^^^.. 

^ Eirst results indicate that this requires large system sizes to reduce 

i-C the uncertainties in the finite-size scaling Ansatze. Nevertheless, our 

investigation shows that within the achievable accuracy the numerical 
estimates are still compatible with analytic predictions, contrary to 
the recent claim by Bock and Vink. 
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1 Introduction 



Two-dimensional (2D) Euclidean quantum gravity is believed to be an im- 
portant toy model on our way to a realistic quantum theory of gravity. Along 
with the theoretical achievements of string theories 2D Euclidean quantum 
gravity has become a rather well understood subject. In contrast to the clas- 
sical theory, which is dynamically trivial, the quantum theory possesses a 
rather rich structure due to the conformal anomaly. One of the interesting 
aspects is the intrinsic fractal structure of space time which shows up in the 
divergence of the partition function Z[A) with increasing area A, governed 
by the string susceptibility exponent 7str [1]: 



Here A^^ denotes the renormalized cosmological constant. The exponent 7str, 
which depends on the genus g of the surface, has been calculated first by 
Knizhnik, Polyakov, and Zamalodchikov (KPZ), using conformal field theory 
methods [2] to be 



Because 7str = 2 is the classically expected result, quantum effects can only be 
seen for non-toroidal topologies [g ^ 1). Eor spherical topologies [g = 0) the 
prediction 7str = — 1/2 agrees with matrix model methods [3]. Moreover, if 
one couples spin matter to gravity, conformal field theory and matrix models 
both predict a significant change in the critical exponents of the associated 
spin phase transition when compared to fiat space exponents [2, 4]. These 
remarkable agreements strengthened the belief in the existence of universal 
features of 2D quantum gravity. 

The matrix model approach gave rise to a novel numerical scheme, called 
the dynamically triangulated random surface (DTRS) approach to quantum 
gravity [5]. The DTRS method could reproduce the prediction for 7str [1] as 
well as the critical spin exponents [6]. 

Eor the alternative numerical approach using Regge calculus [7, 8] the 
situation is far away from being clear. One major conceptual problem is 
associated with the path integral measure of quantum gravity (for a good 
review see Ref. [9]). There have been almost as many different proposals for 
the integration measure as there are different approaches to quantum grav- 
ity. Within the gauge theory approach, depending on one's own preference 
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in the correct gauge group of gravity, such as diffeomorphism, local Lorentz, 
Poincare or conformal group, one obtains different results for the integration 
measure. The same situation holds true if one uses the Hamiltonian formal- 
ism or tries to approach the problem by measure theory [9]. Most measures 
differ by the power of the determinant of the metric (/, which stands in front 
of the integration measure, but also non-local objects of the Faddeev-Popov 
type can appear. 

We would like to stress that this ambiguity resides in the continuum ap- 
proach to quantum gravity, and hence should haunt any quantum theory of 
gravity; it is not a problem confined to Regge calculus. In light of the agree- 
ment of conformal field theory and matrix models, however, one is inclined 
to think that at least in two dimensions there might exist something like a 
unique measure and that both, the KPZ and the matrix model approach, 
capture all the essential features of this measure to fall into the same univer- 
sality class. To support this argument one would like to see Regge calculus 
as a third unrelated candidate for a quantum theory of gravity to repro- 
duce all known results (for a different opinion, however, see Refs. [10, 11]). 
This would enforce the hope that also in higher dimensions the details of the 
numerical approach are not overly important. 

Unfortunately the results of numerical investigations using Regge calculus 
have been disappointing so far. Using the commonly employed dl jl measure 
on a Regge lattice, no change in the phase transition of an Ising model coupled 
to gravity was observed [12, 13], the critical exponents remained in the flat 
space Onsager universality class. Still there is the hope that with a different 
measure or a different spin coupling to gravity one can reproduce the KPZ 
critical exponents [14, 15]. For pure gravity and toroidal topology agreement 
with 7str = 2 was reported [12, 16], but, as noted before, the torus is not a 
good testing ground for (2). For the sphere and the dljl measure. Gross and 
Hamber [12] found weak numerical evidence for 7str = —1/2, while Bock and 
Vink in Ref. [16] clearly did not see the dependence of (2). Soon afterwards, 
however, it was claimed in Ref. [17] that the string susceptibility exponent 
7str strong i?^-gravity was consistent with the theoretical predictions [18] 
for the sphere [g = 0), but not for the bi-torus [g = 2). In light of these 
contradicting results, we found it necessary to reinvestigate this problem 
again. 

In the Regge approach part of the difficulty to estimate 7str is the lack 
of methods to measure 7str directly. All approaches used so far introduce a 



2 



curvature square term and deduce from its expectation value an estimate 
on 7str through a finite-size scahng analysis. The problem with i?^-gravity is 
that there are actually two distinct scaling regimes, each of which defines its 
own string susceptibility. For the model defined by 

Z{A) = I V^{g)e-'-5{ j <fx^ - A), (3) 

with the gravitational action taken as 

Sg = j <Px^{\ + -^R'), (4) 

and T>fj,[g) being the DeWitt measure, the scaling behavior of Z[A) was 
investigated by Kawai and Nakayama [18]. The coupling constant a has 
the dimensions of an area, and therefore sets a length scale of -^/a. The 
dimensionless quantity A := A^a can then be used to distinguish between 
the cases of weak i?^-gravity (A ^ 1) and strong i?^-gravity [A <^ 1). For 
the case A ^ 1 the scaling relations (1) and (2) are recovered, whereas for 
A ^ 1 it was found, that 

Z{A) oc A"'--\-^=/^e-^«^-«^, (5) 

where 

7^ = 2-2(1 -^/), (6) 

Sc = 167r^(l — gY is the classical action, and ^ is some constant. Only for 
the torus [g = 1) the two scaling behaviors and string susceptibilities are 
the same. Note that the dependence of 7^^^. on the gender of the surface is, 
loosely speaking, weaker than that of 7str- For later use we also write down 
the derivative of Z with respect to A: 

- a^R (A > 1), (7) 



dA A 

--a\R +SJA^-li (A«l). (8) 



d\nZ 7^tr 



dA A 

While the partition function Z is not directly accessible in Monte Carlo 
simulations, logarithmic derivatives like (7) and (8) can be estimated by 
measuring appropriate expectation values. 
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The remainder of the paper is organized as follows. In Sec. 2 we first re- 
view the Regge discretization method. Then we discuss the finite-size scaling 
methods used in Ref. [16, 17] to determine 7str and 7^^^., show its shortcom- 
ings, and suggest a new method to measure 7str and 7^^^.. Section 3 deals 
with some details of our Monte Carlo simulations. In Sec. 4 we present our 
numerical results and discuss their interpretation, and in Sec. 5 we conclude 
with a brief summary of the main results and some final remarks. 

2 Discretization Method and Finite-Size Scal- 
ing Ansatze 

2.1 Regge discretization 

Regge's discretization program [7] consists of replacing a given continuum 
manifold by piecewise linear manifolds, whose internal geometry is fiat. This 
procedure works for any space-time dimension and for metrics of arbitrary 
signature. Here we restrict ourselves to the simplest case of two dimensions 
and Euclidean signature. 

In two dimensions this procedure is most easily visualized by choosing 
a triangulation of the surface under consideration, where each triangle then 
represents a piecewise linear manifold. The net of triangles is itself a two 
geometry, with singular (non-differentiable) points located at the vertices of 
the net, where several triangles meet. A vector that is linearly transported 
around these vertices experiences in the presence of curvature a rotation by 
the deficit angle 5i = 2tt — Y^tDi ^ii^)} where 6i{t) is the dihedral angle at 
vertex i. For the 9jrG9j 9jSS1 gnment we used the barycentric decomposition, 
where Ai = Y^f^i At/ 3 denotes the barycentric area with At being the area of 
a triangle t. We can then identify the following continuum quantities with 
their discrete counterparts, namely 




(9) 




(10) 
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Our notation is identical to that used in Ref. [13]. In two dimensions the 



Einstein-Hilbert action / cPx^J g{x)R{x) is by the Gauss-Bonnet theorem a 
topological invariant, which makes such a theory classically trivial, there are 
no equations of motion. Regge [7] gave a beautiful proof of this theorem in 
terms of the deficit angle. The sum over the deficit angles in two dimensions 
is proportional to the Euler characteristic, namely ^j ^i = 47r(l — g). We 
therefore will not consider the corresponding term in the gravitational action. 
If we keep the area A fixed to its initial value, then, classically, dynamics 
can only arise from the i?^-interaction term. Such a term was used in two 
and higher dimensional studies to cure the unboundedness problem of the 
gravitational action [19]. 

For triangulated surfaces the Euler relation reads as 



where A^o,^i, and N2 denote the number of sites, links and triangles, re- 
spectively. Eor triangulations without boundary we also know that a link is 
shared by two triangles, resulting in the relation = N2I2. Erom these 

two relations one can derive two more, namely A^o — 2(1 — g) = N2I2 and 
A^o — 2(1 — g) = which will become useful later. 

Eor each triangle there is a one-to-one correspondence between the square 
of the link lengths and the components of the metric. Denoting by g^u{i) the 
components of the metric tensor for the i^^ triangle, and by qi^^^i^^^qi^i^n, 
and Qi^i+i, the square of its three edge lengths, one can derive the following 
relation g^ipii) = \ + — In classical Regge calculus one 

starts with the action principle and derives the equations of motion, one for 
each link. The link lengths have to be adjusted to satisfy those equations in 
order to be a classical solution. The connectivity of the edges, in simplicial 
topology called the incidence matrix, is fixed from the beginning through the 
simplicial decomposition of the manifold under consideration. 

In quantum Regge calculus the technical aspects are similar, although 
the philosophy is quite different. Here we want to evaluate the functional 
integral in Eq. (3) by Monte Carlo (MC) methods. In principle, the integral 
has to be extended over all metrics on all possible topologies, but, as usual, 
we restrict ourselves to a specific topology, here the sphere. The integral 
over the metric is replaced by an integral over the square of the link lengths. 
An important ingredient in the functional-integral method is the appropriate 
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measure, which, as already explained in the Introduction, is not even known 
in the continuum. The most popular measure is DeWitt's supermetric [20], 
a distance functional on the space of metrics. It was used by Polyakov in his 
famous string solution [21]. Because in 2D the measure is the primary source 
of the non-trivial dynamical content of the theory, its correct transcription 
might be the key point in reproducing the KPZ results. Nevertheless if 
the discretized DeWitt measure is still a local one, then one might argue 
on the basis of universality that any other local measure will do as well. 
In the present study we only report simulations with the most commonly 
used "computer" measure T>ji{l) = Yl(ij) T^Fih])-, in order to facilitate a 
comparison of our data with the results of Ref. [16]. Here F{1) is a function 
which takes the value one, if the links /ij, forming one triangle, fulfill the 
triangle inequalities, and zero otherwise. For small link lengths we introduced 
a scale invariant cut-off to avoid round-off problems. We will take up the 
question of the measure dependence of 7str (7str) ^ forthcoming publication 



Collecting the transcriptions from the continuum to the simplicial Regge 
approach, our lattice analogue to Eq. (3) is therefore given by 



where we abbreviated R^^ := jAi. 

2.2 Finite-size scaling 

The method to extract 7str through the finite-size scaling (FSS) properties of 
the expectation value of an added i?^-interaction term was first discussed by 
Gross and Hamber [12], and later improved in Ref. [16]. Here we would like 
to point out that there is still room for ambiguities in the FSS assumptions. 
A very simple derivation of the scaling behavior comes from a rescaling argu- 
ment. We will not discuss in this paper any measure ambiguities, and restrict 
ourselves to the scale invariant computer measure o?///, as used in our (and 
most previous) simulations. Other measures could be treated without any 
problem. The simple derivation presented below has the additional advan- 
tage over the one given in Ref. [16], that it neither invokes Ward identities 
nor any artificial weight functions. 



[22]. 




(13) 
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Consider the partition function (13). A rescaling of the hnk lengths / to 
the dimensionless link lengths /' of the form / — y 1'\/A yields Ai — V 
AY,i A'-^ and therefore 



From this expression one can easily compute the derivative 

dlnZ ,1 



aA + ^(i?'-l), (15) 



dA A 
where 

R':=a{j:Rj) = a{j:5f/A,), (16) 

i i 

and by inspecting (14) it is also easy to see, that = i?^(A, A^i), meaning 
that the expectation value depends only on A^i and the dimensionless param- 
eter A. For large A^i one expands the finite part of R^ into a power series, 
whose first three terms read as 

= ... + boA + W + b^/A + .... (17) 

If this is inserted in (15), then a comparison with (7), (8) gives bo = —{Xr — 
A)a, 6i = 7str - 2 for A > 1, and bo = - {Xr - A)a, bi = 7^^^. - 2, 62 = Sc 
for A < 1. 

The scaling Ansatz of Refs. [16, 17] is to consider first a power-series 
expansion of R^[A^Ni) in A^i. They actually use N2 instead of A^i, but this 
is trivial because for any compact triangulation we have the fundamental 
relation 3A^2 = 2A^i. Furthermore, their expansion of R^ is not done at 
fixed A, but at a fixed discretization scale set by the average triangle area 
ao := A/N2 through the dimensionless parameter cio := ao/a: 

R\ao, N2) = N^coiao) + Ci(ao) + C2(ao)/A^2 + • • • • (18) 

The coefficients Ci[ao) are thus defined in the thermodynamic (infinite area) 
limit. This expansion has to be justified by the simulation results, because 
it is not based on any ah initio calculations. A similar good educated guess 
would have been to to consider for instance an expansion in a linear length 
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parameter ^/W^. In a second step the coefficients q are expanded in Refs. [16, 
17] into a power series in ciq as 

Co = 4°) + ao4') + ..., (19) 

Cl = cS°^ + aocS'^ + . . . , (20) 

C2 = (cf + aoc^') + . . .)/ao, (21) 

and the continuum hmit is taken by sending the discretization scale to zero, 
ao — > 0. A comparison with (17) then yields ho = Cq^\ bi = c^°\ and 62 = C2°^ 
Note, that in order to make contact with the continuum result of Eq. (8), 
C2 needs to start with a divergent term I/ciq. Only in the combined limit 
N2 — y oo,ao — y 0, this makes sense. But because ciq is fixed, and not A, 
effectively there is no control of the crossover from strong i?^-gravity scaling 
behavior (A <^ 1) to the weak i?^-gravity scaling behavior (A ^ 1). If one 
takes first the thermodynamic limit in (18) then one always obtains the values 
of the coefficients Ci in the limit A = A^20o ^ 1, hence for weak gravity. 
This means, c|°^ = 7str — 2, and C2 — > 0, but one has to be careful to make 
the system size always large enough to reach this limit. If one truncates the 
fit at some suitable value of N2 to explore the region where A ^ 1, then the 
continuum limit can only be taken at finite N2. Because the results of Ref. [17] 
were obtained in for very small N2 finite-size effects can become important. 
Even worse, because the coefficients cp') in the expansion (20) are constants, 
how can c|°^ + 2 change from 7str to 7^^^., as is claimed in Ref. [17]? These 
subtleties, which have not been previously addressed, make it appear very 
unlikely that one can unambiguously extract 7^^^. in this way. 

We therefore suggest an alternative approach, where we look at the ESS 
behavior at a constant value of A. This is very much in the spirit of the 
investigations of Ref. [18]. Expanding i?^(A, A^2) at constant A we obtain 

i?2(i, N2) = N2do{A) + (/i(i) + d2{A)/N2 + .... (22) 

Eor large A we have to increase N2 to insure that the following terms remain 
small. Also here other ESS Ansatze are possible. The next step is to expand 
the coefficients di as a power series in A. This time the coefficient di carries all 
the necessary information to extract the string susceptibilities. A comparison 
with (17) yields 

(/i(i)= 60A + (7,tr - 2) + 0(l/i) fori>l, (23) 

diiA)= SjA + iY^,^-2) + boA + OiA^) for 1 « 1. (24) 
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If we plot 0?! versus A we thus expect to see a linear behavior for very large 
A and a divergent behavior for small A, governed by the classical action 
Sc = 167r^(l — gy. The difference between our method and that of Ref. 
[16] appears as a subtle interchange of thermodynamic and continuum limit. 
We first take the continuum limit (A^2 — ^ oo) for fixed A, and then the 
thermodynamic limit, whereas in (18) first the thermodynamic limit is taken 
for fixed ciq = A/A^2, and then the continuum limit is performed. 

The appearance of the classical action is not hard to understand. Actu- 
ally, for any regular triangulation with coordination number q of arbitrary 
topology we find from the Euler relation 5i = 47r(l — g)/q^ and Ai = A/q^ 
therefore ai?^ = aY^5f/A, = 167r^(l — gY / A. For small A this will be the 
dominant term. 

Let us assume that we have a distribution of the square of the link 

lengths that give the classical action Sc = For small A one should 

be able to make a stationary (semi-classical) approximation of the action 
S{q) around S'(g^°^) in the variables g, which should become small in the 
continuum limit due to the constancy of A. Our results suggest that it looks 
like 

S{q)^Sc + NoF{A)+ Y: ^'V°^)(%--g!?)(g./-gi?) + ..., (25) 

<ij><kl> 

where the primes denote differentiation with respect to the g, and S''(g^°^) = 
0. F denotes the first correction to Sc independent of g, which is responsible 
for the divergent behavior of proportional to A^o- This form of (25) would 
heuristically explain how do and the classical action part in di arises. To 
zeroth order always the classical part dominates, and it is only in the limit 
of large A that quantum fiuctuations become more important. If one could 
compute the second derivative part of S{q)^ then one could probably also 
calculate 7^^^.. 

Notice that if one uses A^o instead of N2 in Fqs. (18) and (22), one will 
get in general different coefficients q and c/i, respectively. An expansion of 

i?2(ao, A^o) = Noel + cl + c^No + . . . (26) 

instead of (18) would give for spherical topologies Cq = 2co, c[ = Ci — 4co, 
and C2 = C2/2. A further extrapolation of the coefficients q to cio — > would 
therefore yield the same string susceptibility only if Cq vanishes in this limit. 
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This is not supported by our data, where we see that Cq actually increases 
as ciq decreases. For the case of Eq. (22) the situation is similar. The use of 
A^o instead of N2 results in c/q = 2c/o, = di — 4c/o, and c/2 = 0^2/2. This 
results in a decrease of the string susceptibilities by Ado. We stress that we 
believe that, a priori, should be considered as a function of A^i, and that 
it should therefore be expanded in A^i. Because the use of N2 is just a simple 
scale change, this would be equally valid. Still it is troublesome that even 
in the continuum limit, the difference between the expansions using N2 and 
A^o remains. From the Fuler formula one can see how the difference depends 
for general topologies on the gender of the surface. The only changes are 
cl = Ci — 4(1 — g)co and dl = di — 4(1 — g)do. Again only for the torus 
(g = 1) this ambiguity is absent. 

We attribute above ambiguity in the determination of 7^^^. (7str) to the 
fact that 7gjj. (7str) comes from a sub-leading effect, whereas the leading effect 
Z oc 6""^=/^ can be determined with both, A^o and N2. For Fq. (18) this can 
be seen by simply inserting the definitions of the coefficients. One then finds 
^(0)* _ 1^(0) _ Pqj, (^22) one needs to assume that do{A) has no l/A 

dependence for A <^ 1. 

3 The Simulation 

All MC simulations were performed on spherical topologies, which were re- 
alized as triangulated surfaces of a three-dimensional cube, as discussed in 
Ref. [16]. For this choice of discretization six vertices have coordination 
number four, whereas the rest has coordination number six, see Fig. 1. 
The number of vertices A^o is related to the edge length L of the cube by 
A^o = 6(iv — 1)^ -\- 2. Usually, the size of the lattices varies from L = 7 up 
to 55, corresponding to 218 up to 17498 lattice sites, or 648 up to 52488 
link degrees of freedom. To update the links we used a standard multi-hit 
Metropolis update with a hit rate ranging from 1 ... 3. In addition to the 
usual Metropolis procedure a change in link length is only accepted, if the 
links of a triangle fulfill the triangle inequalities. 

The area A was kept fixed at its initial value A = Y^i Ai = N2/2 during 
the update to simulate the delta function in Fq. (13). In order to achieve 
this, we would need in principle to rescale all links during each link update, 
amounting in a non-local procedure. However, due to the scaling properties 
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of the partition function, this can be absorbed in a simple scale factor in 
front of the term. To avoid round-off errors we explicitly performed a 
rescaling after every full lattice sweep. Notice, that technically our simulation 
procedure is different from the methods employed in Refs. [16, 17]. 

The first set of simulations was designed to test the extrapolations to the 
continuum limit which in Ref. [16] were based on very small values of I/ciq in 
the range of 0.5 — 5. Already exploratory runs [23] indicated that this range 
is not appropriate to see the true asymptotic behavior. In the present study 
we therefore used much larger values of I/ciq = 0.2, 2.5, 10, 20, 40, 80, 160, 
320, 640, and 1280. 

The second set of simulations consisted of runs at constant A at the 
values of i = 9126/a with a = 1.25,2.5,5,10,20,40,80,160,320, and 640, 
which covers roughly the range of A = 14 — 7300. 

For each run we recorded about 10000 measurements of the curvature 
square = Y^iSf/Ai. To reduce temporal correlations, we took measure- 
ments every second MC sweep through the entire lattice. The statistical 
errors were computed using standard jack-knife errors on the basis of 20 
blocks. The integrated autocorrelation time tji2 of R^ was usually in the 
range of 5 - 10. 

4 Results 

4.1 Scaling at fixed 

We first take a look on the raw simulation data, and plot /N2 versus I/N2. 
The set of data points for the four largest values of ciq can be inspected in 
Fig. 2. Because we chose ciq = l/2a this is equivalent to the simulations 
with the four smallest values of a. The curves look straight to the eye, and 
the scaling Ansatz (18) seems to works well even without the C2 coefficient, 
indicating that N2 is large enough so that we are in the weak i?^-gravity 
regime. 

A closer look on the curves for I/ciq = 10 (l/oo = 20), see Fig. 3, shows 
apparently two scaling regions, divided approximately by a line through 
1/A^2 ~ 0.005(0.003). We interpret this region as the crossover region from 
A ^ 1 to A ^ 1. Because A was chosen to be A^2/2a, we decrease A ei- 
ther by decreasing A^2 or by increasing a. This means we always start out 
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on small lattices in the strong i?^-gravity regime, and end up on sufficiently 
large lattices always in the weak i?^-gravity regime. It is therefore hard 
to imagine, that one can fit the whole range of data points with the same 
truncated Ansatz (18) without taking into account contributions from other 
coefficients which arise from interchanging the A limits. We know the scaling 
behavior only for the two limiting cases of A, but nothing in between these 
two limits. By increasing N2 we would again be able to extract Ci with a 
linear Ansatz, yielding 7str, but already at the rather coarse discretization 
scale of ciq = 0.05 — 0.1, the studied system sizes turned out to be too small 
to produce a reliable estimate for Ci. 

For the lower values of cio, the data points of the smaller lattices shown 
in Fig. 2 begin to show a clear deviation from the straight line behavior. In 
Fig. 4 we show the raw data for W j N2 for the four smallest values of ciq. On 
the basis of Ansatz (18) this means, that the C2 coefficient, originating from 
the exponential damping factor of the classical action Sc for the case A ^ 1, 
wins more and more. To make this effect more visible and to get more data 
points for small A, we performed four additional simulations for the four 
largest values of a on the lattice sizes as small as A^o = 26 — 152 (L = 3 — 6). 
We now mimic the procedure of Ref. [17] to extract an estimate for 7^^^. by 
staying at a rough discretization scale to achieve a sufficiently small A. 

We fitted all data points using Ansatz (18), including a hnear and a 
quadratic term in 1/A^2, see Table 1. In Ref. [16] it is shown, that in the 
hmit of ciq — > 0, Ci should approach 7str — 2. Here, because we truncated 
our fit at some A^2, we approach this limit by increasing a, and therefore 
decreasing A, so that effectively Ci should approach 7^^^. — 2. But, as already 
mentioned before, this procedure is really ill defined, because it is simply 
impossible to take the continuum limit of ciq — > 0, and staying at the same 
time at A ^ 1. We are therefore not convinced, that the coefficient Ci is 
simply related to 7^^^.. 

Nevertheless, by recalling Fqs. (17)-(21) the value of C2ao should approach 
167r^. As can be inspected in Table 1, C2ao indeed approaches 167r^, but the 
total of fit is unacceptably high. The steep increase in the total 

of the four simulations with highest a of course originates in the higher 
number of data points on the smaller lattice sizes. For those simulations the 
tail of the power law decay, namely the data points on the large lattices, 
produce the increase in x^ ■ would expect then that removing the data 

of the larger lattices would lead to an improvement, because they belong to 
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the regime where A is large and the C2 coefficient is not so important. We 
therefore successively discarded the larger system sizes until we obtained a fit 
with an acceptable quality. The final values of the fit parameters along with 
the remaining number of degrees of freedom (dof) can be found in Table 2. 
Indeed we see that we can stabilize the fits much better, and that the value 
for the classical action Sc comes out very well. Only the simulations with 
largest coupling a show the first signs of numerical problems we experience 
with large values of a. 

Now we have the paradoxical situation that one needs very small lattice 
sizes or very large values of a to expect good results. On very small system 
sizes, however, one needs to worry about finite-size effects, which can lead to 
rather large systematic errors. As a final estimate of this analysis we take the 
average of the four MC estimates, but due to the aforementioned criticism 
of this method, we enlarge the error, leading to 7^^^. = c^°^ + 2 = 0.1(3) as a 
conservative estimate. From (6) we see that this value is still consistent with 
the prediction 7^^^. = 0. 

In Ref. [17] a similar value was found for the sphere, but the value of 7^^^. 
for the bi-torus [g = 2) did not agree with the theoretical expectations. We 
suggest two possible explanations for this. The first is that the presence of 
large FSS corrections is responsible for the "failure" of the bi-torus analysis, 
because the results were obtained on relatively small lattice sizes. The other 
is that actually the methods presented in Ref. [17] cannot predict anything 
about 7gjj., because inherently the continuum limit in ciq cannot be taken. 
Therefore the coincidence of our results with the theoretical value of 7^^^. for 
the sphere might be purely accidental. 

Still we investigated our data a bit more closely. Another way of improv- 
ing the quality of the fits is to try out other scaling Ansatze or to include 
more correction terms to (18). We therefore used (18) augmented by a term 
C3/NI: 

R^cio, N2)/N2 = co(ao) + ci(ao)/A^2 + C2{do)/N^ + C3{do)/Nl (27) 
We also tried an Ansatz of i?^(ao, A^2) in terms of ^/N2- 

R^{do, N2)/N2 = co(ao) + ci(ao)/A^2 + C3/2(ao)/A^2^' + C2(ao)/A^2' + • • • • (28) 

The results for the two generalized fits can be inspected in Tables 3 and 
4. We find that both Ansatze seem to work equally good or bad, so that 
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on the basis of these results alone one cannot draw any conclusion on their 
validity. Actually, the inclusion of more correction terms does not improve 
the simulations with large and small values of a, only the crossover region 
seems to get improved. Finally one can fit the values of the coefficient Ci 
obtained according to the Ansatze (18), (28) and (27) versus 1/ao, see Fig. 5, 
to make the extrapolation ciq — y and to see how the observations of 
Ref. [16] comply with our findings. Fig. 5 was composed using the values 
from Tables 1, 3, and 4. For small values of ciq all three curves seem to come 
closer together, and are approximately around the theoretical value of 7^^^.. 
For large values of ciq one sees no way how to extract 7str- If one extrapolates 
to large values of cio, then this means A ^ 1, but in the same limit the 
discretization scale becomes very large, and we can not be sure if we still 
reach the continuum limit. 

To conclude this subsection it should be stressed that at the couplings 
where we can compare with with Ref. [16], the raw data for do agree 
within error bars, so that differences in the final results cannot be blamed 
on using different simulation techniques. It is rather our much larger lattice 
sizes and the considerably increased range of ciq which reveals the potential 
problems with the approach of Refs. [16, 17]. 

4.2 Scaling at fixed A 

The raw data of our simulations at fixed A is shown in Fig. 6. One first 
notes that all curves of /N2 versus 1/A^2 for different fixed values of A are 
significantly curved. Straight line behavior is visible only asymptotically for 
large values of N2. 

Because of the non-linear FSS behavior observed in Fig. 6, the task of 
extracting the coefficient di of Ansatz (22) proves to be a difficult one. In 
our attempt to fit all available data points we therefore used besides the FSS 
Ansatz (22) also a three-parameter fit of the form 

R'{A, N2)/N2 = do{A) + d,{A)/N2 + d^i2{A)lNl'^ + . . . , (29) 

which yielded a better fit, with a that was mostly only half as high as 
for the Ansatz (22). Still neither of those fits yielded a satisfactory total 
■ This can mean that either both Ansatze plainly do not work, or that 
FSS corrections are still so large that one would need even more correction 
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terms, which, in hght of the few available data points, is not applicable. We 
therefore tried to discard the data on the smaller lattices until we obtained a 
fit with a reasonable ■ The values of the fit parameters together with the 
total s-iid the remaining number of degrees of freedom (dof) can be found 
in Tables 5 and 6. As a general trend one observes that the acceptable fit 
range increases with decreasing values of A, with the exception of the two 
simulation with smallest A. The first part of the observation can be inferred 
from Fig. 6, because the curves show more linear behavior the smaller A 
gets. The simulations for the two smallest A do not follow this trend because 
the large values of the coupling a produce numerical problems on the 
larger lattices. Some configurations seem to freeze in local minima and our 
algorithm does not seem to be able to relax those minima sufficiently fast. 

The resulting plot of both values for di versus A is shown in Fig. 7. 
We observe that qualitatively both curves fulfill the theoretical expectations, 
namely they show a divergence at small A and a fiattening slope at large A. 

To extract 7str we need to look at large values of A. Because the divergent 
term goes like IGtt^/A this means large compared with IGtt^. Looking at the 
five highest values of A we note that no clear linear slope can be observed, so 
that with the present data it is impossible to obtain a reliable estimate for 
7str- From a crude linear fit it appears that 7str is too negative, which goes 
just in the opposite direction of what was claimed in Ref. [16]; however this 
observation has to be taken with great care. As far as the extraction of 7str 
is concerned we must thus conclude that either A is not yet large enough to 
see the linear behavior or the presence of large FSS corrections in the first 
fit have led us to underestimate the errors in our data. The first obstacle 
could be overcome by increasing A, whereas the second would require even 
larger system sizes. Unfortunately it turns out that we need to follow both 
suggestions simultaneously to improve our data. We have seen in Fig. 6 
that with increasing A the curves rapidly deviate away from linear behavior 
already for quite large system sizes. Therefore increasing A requires also 
increasing the system sizes to obtain a comparable accuracy, which in turn 
requires more and more computing efforts which are beyond the scope of the 
present study. 

The situation for 7^^^. appears somewhat better. Performing a three- 
parameter fit of the form (24) to the four data points with A < IGtt^, we 
obtain Sc = 399(46) and 7^^^. = —2.2(1.9), if di is taken from the fits to the 
Ansatz (22), and Sc = 608(77) and 7^^^. = —6.6(3.1), with di from Ansatz 
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(29). Already the deviations of the two fit results from each other show that 
one should not trust these numbers. At least for Sc we know that a value of 
IGtt^ 157.91 should emerge. This shows that in principle we should include 
more FSS corrections to stabilize the fits, apart from the obvious need for 
more data points. 

The next thing we tried is to discard in the first fits at fixed A the data 
points for small systems to reach a range where a linear fit of the form 

R^{A, N2)/N2 = do{A) + (/i(i)/A^2 (30) 

is sufficiently accurate. The resulting fits for the six smallest values of A are 
shown in Fig. 8 and the fit data can be found in Table 7. For sufficiently large 
N2 we see a linear behavior, only the two simulation sets with smallest A are 
somewhat scattered around the fit line. As already mentioned we attribute 
this to numerical problems in thermalizing the configurations with extremely 
large coupling, which is presumably not refiected by the computed jack- 
knife errors. The fit over the resulting six data for di according to Fq. (24) 
is shown in Fig. 9a and yields Sc = 199(65) and 7^^^. = 5.1(3.0), with a 
reasonable total chi-square of 3.5. By looking at the graph, however, 

one gets the impression as if the four data points with lowest A lay more 
on a scaling curve like Fq. (24), suggesting that the values of A > 158 are 
already too large to show the asymptotic behavior. Indeed, these values may 
already belong to the crossover regime to the weak scaling. In Fig. 9b 
we therefore fitted just the first four data points, resulting in Sc = 187(104) 
and 7gjj. = 5.1(6.6), with = 0.63. To see what difference it would make to 
constrain Sc to its theoretical value of 167r^ we also included a two-parameter 
fit of the form di = 1671^ A + 7^^^ -2 + boA, that resulted in 7^^^ = 7(7). To 
the eye both curves are almost identical. Within the relatively large error 
bars the values of both fits are compatible with the theoretical expectations. 
In comparison with the analysis in the preceding subsection the estimates 
appear less accurate, but this time uncontrolled systematic finite-size effects 
are definitely reduced. To enhance the accuracy we would again need much 
more data points, now at small values of A. Here it not necessary to go to 
very large system sizes in the first fit, but one faces the numerical problem of 
simulating at extremely large R^ coupling, where the generation of a reliable 
sample of MC configurations for probing the partition function in (13) proved 
to be very difficult. 
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5 Conclusions 



We tried to measure the string susceptibilities 7str and 7^^^. using two FSS 
Ansatze in . Although the approach of Ref. [16] is in principle applicable 
to determine 7str correctly, our results for the dl/l measure show that it fails in 
practice because the system sizes needed in order to reach the thermodynamic 
(large area) limit N2 — > 00 on reasonable small discretization scales Oq, are 
too large. Effectively this can be seen in a crossover behavior from weak 
to strong scaling, which is mainly due to the fact, that already in the 
first truncated Ansatz of Eq. (18) one mixes data with small and large A by 
varying N2. Only if N2 could be made sufficiently large, one would be able 
to stay always in the large A regime. 

The method of Ref. [17] to estimate 7^^^. with the same Ansatz can be 
performed practically, but here we encounter severe conceptual problems. 
The analysis was done employing only very small lattice sizes, and to reach 
a sufficiently low A one has to discard successively the larger lattice systems 
and to extrapolate to very small system sizes. In this way one will neces- 
sarily experience large finite-size effects. Our estimates for 7^^^., obtained for 
smaller discretization scales and larger lattices, agree with the finding of [17], 
showing, that finite-size effects seem to be unimportant for the sphere. Nev- 
ertheless the continuum limit — > can only be performed for some finite 
value of N2. The coincidence of our value for c|°^ -|- 2 with 7^^^. might there- 
fore be purely accidental and should not be considered gnificant test 
of Regge calculus. In light of all these conceptual problems the "failure" of 
the bi-torus analysis in Ref. [17] should not be regarded as a serious problem 
for Regge calculus, but rather as a problem with the method itself. 

Alternatively, if one works at a well controlled A as in (22), one should in 
principle be able to explore both scaling limits and to predict both 7str and 
7str- have experienced, however, large ESS corrections in the weak R^ 
regime, and numerical problems in the strong R^ regime, which make it very 
difficult to extract the coefficient di with high accuracy. Unfortunately, these 
are the regions in which a precise knowledge of di[A) is needed to extract 
the asymptotic behavior of Z[A) and the associated string susceptibility 
exponents. Eor this reason we were not yet capable to measure either 7str or 
7gjj. with high enough precision to decide about a "failure" or "non-failure" 
of Regge quantum gravity. The results we obtained so far are still consistent 
with the theoretical prediction of Ref. [18] for both 7str and 7^^^., although we 
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must admit, that the present data does not provide compeUing evidence for 
either direction. Conceptually, however, already these exploratory results are 
encouraging and we believe that with higher computational effort one can 
arrive at a conclusive statement. A particularly challenging test would be 
the topology of the bi-torus or surfaces of higher gender where the methods 
of Ref. [16, 17] seem to fail also for 7^^^.. 

Still, in the long term, it would of course be desirable to develop more 
direct approaches to measure 7str, as has been done for the DTRS method 
[1]. There has been one alternative method to test the scaling behavior of 
Regge gravity [24], where the universal loop length distribution is used. On 
the basis of a numerical study it is shown, that the dljl measure gives the 
theoretically predicted distribution at least for the baby-loop distribution, 
but only for very large system sizes. Unfortunately it is not yet known, if 
the scaling of the loop length distribution is directly related to the string 
susceptibility 7str- However, for the uniform dl measure, no agreement was 
found. This naturally raises again the question of the existence of a correct 
measure for quantum Regge calculus [15, 23, 25]. Other, more physically 
motivated measure choices are presently under investigation. We conclude 
with the main result of our study, that even for the dljl measure, a failure 
of Regge calculus to produce the theoretical predictions about the string 
susceptibilities has not yet been shown. 
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a 


Co 


Cl 


€2(10 




0.10 


0.1093(1) 


1.7(4) 


-871(766) 


0.8 


1.25 


0.2106(1) 


3.1(2) 


-26(26) 


0.5 


5.00 


0.2332(1) 


4.8(2) 


27(11) 


9.0 


10.00 


0.2391(1) 


4.6(3) 


94(7) 


52.2 


20.00 


0.2430(1) 


5.9(5) 


104(6) 


12.1 


40.00 


0.2455(1) 


7.9(7) 


119(4) 


20.0 


80.00 


0.2485(1) 


-0.57(9) 


157.95(5) 


239 


160.00 


0.2495(1) 


-1.00(5) 


158.03(1) 


419 


320.00 


0.2509(1) 


-1.18(9) 


158.00(2) 


122 


640.00 


0.2525(1) 


-1.74(9) 


157.99(1) 


152 



Table 1: Fit results for the Ansatz ^V^2 = Co + C1/N2 + C2/iV|. We fitted 
over all available data points. For the four largest values of a we added 4 
runs at system sizes of A^2 = 48 — 300 (L = 3 — 6). As a measure for the 
goodness-of-fit we included the total of the fit. For small A we expect 
02(10 = 16^2 ^ 157.91. 
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dof 


a 


Co 


Cl 


€2(10 


x' 


3 


80.00 


0.2530(4) 


-2.0(2) 


158.42(6) 


11.4 


5 


160.00 


0.2520(2) 


-1.7(6) 


158.14(2) 


16.4 


6 


320.00 


0.2522(2) 


-1.6(1) 


158.04(2) 


17.0 


6 


640.00 


0.2535(2) 


-2.1(1) 


158.00(1) 


75.3 



Table 2: Fit results for the Ansatz R'^/N2 = Co + Cil N2 + C2I N^. We discarded 
successively the largest lattices until we obtained a reasonable total of the 
fit. The number of degrees of freedom for the fit is denoted by dof . For small 
A we expect C2ao = 167r^ 157.91. 



a 


Co 


Cl 


02(10 


C3 


x' 


0.10 


0.10929(8) 


1.7(8) 


-793(5519) 


-5146(356559) 


0.8 


1.25 


0.21063(2) 


3.0(3) 


-6(170) 


-16999(143587) 


0.5 


5.00 


0.23309(4) 


5.6(5) 


-134(68) 


541437(225053) 


3.2 


10.00 


0.23877(5) 


8.6(7) 


-207(43) 


1913421(269835) 


1.9 


20.00 


0.24275(9) 


9(1) 


-28(40) 


1772617(519739) 


0.4 


40.00 


0.24519(11) 


13(2) 


6(26) 


2951788(674917) 


0.9 


80.00 


0.24826(8) 


0.6(2) 


156.4(3) 


10324(1323) 


178 


160.00 


0.24944(5) 


-0.4(1) 


157.7(7) 


4693(758) 


380 


320.00 


0.25091(7) 


-1.1(2) 


157.96(5) 


978(1254) 


122 


640.00 


0.25253(7) 


-2.3(2) 


158.09(3) 


-5396(1345) 


136 



Table 3: Fit results for the Ansatz i?V^2 = Co + Ci/A^2 + C2/A^| + Cs/N^. We 
fitted over all available data points. As a measure for the goodness-of-fit we 
included the total of the fit. For small A we expect C2ao = 167r^ 157.91. 
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a 


Co 


Cl 


C3/2 


€2(10 


x' 


0.10 


0.1093(2) 


1.5(2.0) 


6(105) 


-1288(6981) 


0.81 


1.25 


0.2106(1) 


2.9(7) 


10(41) 


-80(225) 


0.44 


5.00 


0.2330(1) 


7.6(1.2) 


-163(63) 


248(86) 


2.19 


10.00 


0.2386(1) 


14(2) 


-515(77) 


428(50) 


6.65 


20.00 


0.2426(2) 


15(3) 


-519(150) 


279(51) 


0.02 


40.00 


0.2449(2) 


22(4) 


-816(185) 


255(31) 


0.44 


80.00 


0.2480(1) 


4.1(5) 


-98(9) 


160.9(3) 


118 


160.00 


0.2493(1) 


1.4(3) 


-50(6) 


158.8(1) 


337 


320.00 


0.2509(1) 


-0.2(5) 


-20(9) 


158.16(7) 


117 


640.00 


0.2525(1) 


-2.5(5) 


16(10) 


157.92(4) 


149 



Table 4: Fit results for the Ansatz ^V^2 = Cq + C1/N2 + Cg/a/iVa''^ + C2/NI 
We fitted over all available data points. As a measure for the goodness-of-fit 
we included the total of the fit. For small A we expect C2ao = 167r^ 
157.91. 



dof 


A 


do 


0?! 


c/2 


x' 


6 


14.26 


0.2514(2) 


21(2) 


-29143(4579) 


20 


6 


28.52 


0.2508(2) 


4(2) 


-13173(3462) 


28 


8 


57.04 


0.24971(5) 


-8.1(4) 


2017(200) 


6.4 


7 


114.08 


0.24904(6) 


-22.7(7) 


10171(604) 


15 


5 


228.15 


0.24830(7) 


-47(2) 


51059(6436) 


7.6 


5 


456.30 


0.2473(1) 


-80(1) 


93358(2356) 


21 


4 


912.60 


0.2465(1) 


-146(3) 


255536(16586) 


14 


4 


1825.20 


0.2450(1) 


-234(3) 


428585(16607) 


17 


3 


3650.40 


0.2434(1) 


-399(5) 


950519(40977) 


15 


3 


7300.80 


0.2402(1) 


-630(4) 


1698794(31367) 


26 



Table 5: Fit results for the Ansatz R^/N2 = do + c/i/A^2 + (^2/A^| with A 
kept constant. We discarded successively the smaller system sizes until the 
fit reached a reasonable total x^. The number of the degrees of freedom of 
the fit is denoted by dof. 
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dof 


A 


do 


0?! 


d3/2 


x' 


6 


14.26 


0.2513(2) 


30(4) 


-991(162) 


22 


6 


28.52 


0.2507(2) 


9(3) 


516(126) 


26 


8 


57.04 


0.24981(6) 


-10.7(7) 


150(15) 


3.5 


7 


114.08 


0.24917(6) 


-28(1) 


502(30) 


8.3 


5 


228.15 


0.24845(9) 


-58(3) 


1452(181) 


6.2 


5 


456.30 


0.2476(1) 


-100(3) 


2649(157) 


15 


4 


912.60 


0.2469(1) 


-184(5) 


5957(384) 


9.8 


4 


1825.20 


0.2457(1) 


-298(5) 


9992(384) 


6.1 


4 


3650.40 


0.2439(2) 


-479(5) 


16865(394) 


21 


3 


7300.80 


0.2419(2) 


-825(7) 


34719(640) 


12 



Table 6: Fit results for the Ansatz ^V^2 = do + di/N2 + c/3/2/A^2''^ with A 
kept constant. We discarded successively the smaller system sizes until the 
fit reached a reasonable total x^. The number of the degrees of freedom of 
the fit is denoted by dof. 



dof 


A 


do 


0?! 


x' 


6 


14.26 


0.2522(2) 


13(3) 


61 


6 


28.52 


0.2508(2) 


6(2) 


38 


5 


57.04 


0.2494(2) 


-5(3) 


4.7 


4 


114.07 


0.2487(2) 


-13(5) 


4.2 


4 


228.15 


0.2483(2) 


-44(3) 


2.7 


4 


456.30 


0.2476(2) 


-82(3) 


2.9 



Table 7: Fit results for the linear Ansatz W IN2 = c/q + o?i/A^2 with A kept 
constant for small values of A. We discarded successively the smaller system 
sizes until the fit reached a reasonable total x^, compare also Fig. 8. The 
number of the degrees of freedom of the fit is denoted by dof . 
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Figure 1: The lattice realization of a spherical topology as the surface of 
three dimensional cube with L = 4^ Nq = 56, and N2 = 108. 
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Figure 2: Scaling of W IN2 versus 1/A^2 for the four lowest values of l/ao- 
Included are the results of the linear fit Ansatz W {ciq, N2) = A^2Co(ao)+Ci(ao)- 



26 



r'/n '-'"'^ 

2 




0.0005 



0.001 0.0015 

1/N 

2 



0.002 



0.0025 



Figure 3: Scaling of W IN2 versus 1/A^2 for 1/ao = 10 and 1/ao = 20. The 
solid and dashed lines serve as a guide to the eye. 
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Figure 4: Scaling of W IN2 versus 1/A^2 for the four largest values of l/ao- 
Note that the smallest system size was lowered to N2 = 48 (L = 3). The 
dashed curve is the result of a fit of the form (18) over all available data 
points. 
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Figure 6: Scaling of W IN2 versus 1/A^2 for all simulations with constant 
A = 9126/a, where A ranges from 14 to 7300. 
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Figure 7: di versus A, where di was extracted according to the two different 
finite-size scaiing Ansatze (22) and (29), showing aii avaiiabie data points. 
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Figure 8: Scaling of W IN2 versus 1/A^2 for simulations with constant A = 
9126/a, where A ranges from 14 to 456. The straight lines denote the linear 
fits according to Ansatz (30). Only the visible data points were used for the 
fits. 
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Figure 9: di versus A, where di was extracted according to Ansatz (30), 
compare Fig. 8. The sohd hues show three-parameter fits of the form Sc/ A-\- 
7gjj. — 2 + boA over the six (a) respectively four (b) visible data points. The 
dashed curve in (b) is a constrained two-parameter fit with Sc = IGtt^. 
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